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ABSTRACT 



' Received / Accepted 

o 

Aims. We study the polarisation properties, magnetic field strength, and synchrotron emission scale-height of Milky-Way-like galaxies 
I in comparison with other spiral galaxies. 

^ , Methods. We use our 3D-emission model of the Milky Way Galaxy for viewing the Milky Way from outside at various inclinations 

as spiral galaxies are observed. We analyse these Milky Way maps with techniques used to obtain the strength of magnetic fields, 
rotation measures (RMs), and scale-heights of synchrotron emission from observations of resolved galaxies and compare the results 
with the Milky Way model parameter We also simulate a large sample of unresolved Milky-Way-like galaxies to study their statistical 
polarisation properties. 

Results. When seen edge-on the synchrotron emission from the Milky Way has an exponential scale-height of about 0.74 kpc, 
■ which is much smaller than the values obtained from previous models. We find that current analysis methods overestimate the scale- 

Ol height of synchrotron emission of galaxies by about 10% at an inclination of 80° and about 40% at an inclination of 70° because 

of contamination from the disk. The observed RMs for face-on galaxies derived from high-frequency polarisation measurements 
^ ■ approximate to the Faraday depths (FDs) when scaled by a factor of two. For edge-on galaxies, the observed RMs are indicative of 

O ^' the orientation of the large-scale magnetic field, but are not well related with the FDs. Assuming energy equipartition between the 

I I magnetic field and particles for the Milky Way results in an average magnetic-field strength, which is about two times larger than the 

O . intrinsic value for a K factor of 100. The number distribution of the integrated polarisation percentages of a large sample of unresolved 

Milky-Way-like galaxies peaks at about 4.2% at 4.8 GHz and at about 0.8% at 1.4 GHz. Integrated polarisation angles rotated by 90° 
C/^ ■ align very well with the position angles of the major axes, implying that unresolved galaxies do not have intrinsic RMs. 

I Conclusions. Simulated maps of the Milky Way Galaxy when viewed from outside at various inclination angles are the basis for a 

comparison with spiral galaxy maps. They are also helpful to check the accuracy of analysing methods of spiral galaxy observations. 
Unresolved Milky-Way-like galaxies are ideal background sources to investigate intergalactic magnetic fields. However, at 1.4 GHz 
they are typically polarised below about 1% and hence a high instrumental purity of future facilities such as the SKA is required to 
I observed them. 

' Key words. Polarisation - Radio continuum: general - ISM: magnetic fields 

cn ■ 

CO 

. 1 . Introduction unresolved galaxies at 1.4 GHz. Simulations with more realistic 

^5 ■ models are needed to achieve a better insight into the problem. 
Understanding the role ofthe magnetic field in spiral galaxies re- Polarisation observations have also been widely used to 

quires to measure its properties and relation to other components study magnetic fields of resolved nearby galaxies (see Beck 2005 

of the interstellar medium. The small-scale magnetic fields of the for a review). RMs are usually obtained from polarisation mea- 

Milky Way Galaxy can be studied in detail. However, properties surements at short wavelengths, such as A3.6 cm and A6 cm, and 

^ of the large-scale field are difficult to recover, since the Solar are believed to indicate the intrinsic galaxy properties. Examples 

. System is located within the disk. are the low-inclination galaxy NGC 6946 (Beck 2007) and the 

Rotation measures (RMs) of polarised background sources "^^1^ ^^ge-on galaxy NGC 253 (Heesen et al. 2009). Here it 

are important probes of the large-scale regular magnetic fields ^as been presumed that galaxies are optically thin for polar- 

in the Milky Way (e.g. Brown et al. 2003, 2007; Taylor et al. f^'^ emission at these wavelengths. However, we have recently 

2009; Van Eck et al. 2011). To obtain a very dense grid of ^0""^ from the Sino-German A6 cm polarisation survey of the 

RMs is the major aim of the Magnetism Key Science Project Galactic plane (Sun et al. 201 1), that the polarisation horizon is 

(Beck & Gaensler 2004) to be cairied out with the future Square ^^out 4 kpc towards the very inner Galaxy at ^6 cm, just a small 

Kilometre Array (SKA) and the various pathfinder telescopes Po^ion of the Galaxy. This motivates us to investigate whether 
such as ASKAP (Johnston et al. 2007) at around 1 .4 GHz. At ^^^^ "^^'"^3' galaxies correctly show the intnnsic RMs of 

this frequency, however, intrinsic polarisation properties of the g^l^^i^s. 

various classes of background sources are poorly understood. The Milky Way Galaxy is a typical spiral galaxy, whose 

Recently, Stil et al. (2009) investigated 23 resolved nearby spiral magnetic fields have been extensively studied by RMs of pul- 

galaxies observed at 4.8 GHz and made a geometric model of the sars (e.g. Han et al. 2006; Noutsos et al. 2008) and extragalactic 

magnetic field to predict the statistical polarisation properties of sources as mentioned above, as well as by total intensity and po- 
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larisation all-sky surveys (e.g. Wolleben et al. 2006; Testori et al. 
2008; Page et al. 2007). By taking into account all relevant ra- 
dio observations available, Sun et al. (2008) and Sun & Reich 
(2010) proposed 3D-emission models of the Galaxy. Although 
many parameters need to be refined, the models are the most 
robust today as they agree with of the observations of RMs of 
extragalactic sources, the all- sky total intensity map at 408 MHz 
(Haslam et al. 1982), the polarisation map at 22.8 GHz (Gold 
et al. 2011), and the WMAP thermal emission template (Gold 
et al. 2009). 

Based on the SD-emission models, we mimic nearby Milky- 
Way-like galaxies by moving the Galaxy away to any distances 
and rotate it by arbitrary values. The large sample of simulated 
Milky-Way-like galaxies allows us to study the statistical prop- 
erties of their integrated polarisation. We also checked whether 
the standard assumption of energy equipartition is relevant to es- 
timate the strength of the average magnetic field, and to what 
extent the inclination of a galaxy aff'ects the determination of its 
scale-height. 

The paper is organised as follows: the simulations of the 
Milky Way Galaxy and its edge-on view are described in Sect. 2, 
RMs, scale-height, and energy equipartition for resolved galax- 
ies are investigated in Sect. 3-5. The integrated polarisation 
properties of a large sample of unresolved Milky-Way-like 
galaxies are studied in Sect. 6. We summarise our results in 
Sect. 7. 



2. Simulations 

2.1. Galactic 3D-emlsslon models 

Details of the Galactic 3D-emission models used in this paper 
were presented by Sun et al. (2008) and Sun & Reich (2010). 
These models are the basis for the current investigations. They 
are able to reproduce the large-scale sky distribution of RMs of 
extragalactic sources as well as total intensity and polarisation 
all- sky surveys available at several frequencies. 

The magnetic field is separated into a large-scale regular field 
and an isotropic random field. The regular field consists of a 
disk and a halo component. The halo field is toroidal with oppo- 
site signs below and above the Galactic plane. The configuration 
of the disk field can be either axi- symmetric plus reversals fol- 
lowing either concentric rings (ASS+RING) or the spiral arms 
(ASS+ARM), or bi-symmetric. Recent observations of numer- 
ous RMs in the Galactic plane by Van Eck et al. (2011) favor 
the ASS+RING field configuration, which will be used for the 
present simulations. The local regular magnetic field strength is 
2 idG and its scale-height is 1 kpc. The strength of the random 
field is 3 yuG elsewhere. 

The 3D-model of the thermal electron density we used is 
NE2001 (Cordes & Lazio 2002). However, it is essential to adopt 
the recently revised scale-height of the diff'use ionised gas by 
Gaensler et al. (2008), otherwise unrealistic parameters result 
as discussed by Sun et al. (2008). We therefore increased the 
scale-height of the thick disk component in NE2001 from 1 kpc 
to 1.83 kpc and changed the midplane thermal electron density 
from 3.4 x 10"^ cm"-^ to 1.4 x 10"^ cm"^. More details were 
discussed by Sun & Reich (2010). For viewing the Milky Way 
Galaxy from outside, we removed all individual features such 
as clumps and voids included in the NE2001 model. The scale- 
height of the cosmic-ray electron density was set to 0.8 kpc, and 
the local density was taken as 6.4 x 10"^ cm"^. We briefly sum- 
marise the main model parameters used for simulations in the 



present paper in Table 1 and refer for details like spatial varia- 
tions to Sun et al. (2008) and Sun & Reich (2010). 



Table 1. Model parameters. 



regular disk magnetic field 


ASS+RING 

zo = I kpc, Rq = 10 kpc 

local field strength: 2/uG 


regular halo magnetic field 


toroidal, asymmetry with plane 
maximum strength: 2/uG 


random magnetic fields 


homogeneous, strength: 3jdG 


cosmic-ray electron density 


Zo = 0.8 kpc, /?o = 8 kpc, y = -3 


thermal electron density 


NE2001 modified 

Zo = 1.83 kpc for thick disk 


volume of the galaxy 


40 kpcx40 kpcxlO kpc 



Notes. Zo' scale-height, Rq: radial scale-length, y: energy spectral index 
which can be related to the spectral index of synchrotron emission as 
a = (y + l)/2 = -1 (/y oc y^, with ly being the intensity). 



In our Milky Way model, we took the observed local syn- 
chrotron emission enhancement into account, which in turn re- 
duced the large extent of the thick disk of previous models (e.g. 
Phillipps et al. 1981; Beuermann et al. 1985) to account for 
the observed emission at high latitudes. These earlier models 
were solely based on the all- sky 408 MHz map and the local 
synchrotron enhancement was not established at that time. The 
problem is that available low-frequency absorption data from 
optically thick HII regions (e.g. Nord et al. 2006; Roger et al. 
1999), summarised by Fleishman & Tokarev (1995), did not 
constrain the 3D-structure of the local enhancement very well. 
However, meanwhile there are more arguments against the exis- 
tence of an extended intensive thick disk of the Milky Way. The 
low percentage polarisation of the WMAP K-band (22.8 GHz) 
map (Gold et al. 2011) requires a dominating random magnetic 
field. However, the regular magnetic field is constrained to a few 
fiG by the RMs of extragalactic sources at high latitudes. This 
in turn means a significant reduction of the cosmic-ray electron 
density to account for the observed high-latitude synchrotron 
emission, which requires to suppress cosmic-ray electron diflfu- 
sion out of the plane. 

The local emission enhancement was modelled as a syn- 
chrotron spheroid with uniform emissivity of 1 kpc radius, cen- 
tered at / = 45°, b = 0° with 560 pc distance from the Solar 
System (Sun & Reich 2010). It is realised in the models by an 
enhanced cosmic-ray electron density, though one can equally 
achieve the same result by increasing the local random magnetic 
fields (Sun et al. 2008). We kept the local enhancement, which 
has no influence on the global intensity distribution of the galax- 
ies, as a mark of the position of the Solar System. 

2.2. Realisation of galaxies from the Milky Way Galaxy 

The simulations were performed with the HAMMURABI code 
(Waelkens et al. 2009). We revised the code to enable simula- 
tions of galaxies in a Cartesian coordinate system besides the 
standard HEALPix system (Gorski et al. 2005). 

To mimic a nearby spiral galaxy, the Galaxy was put at / = /q, 
and b = bo at 3. distance D, and then rotated by Euler angles 
(^1, 02, O3). The process is explained in Fig. 1, where a galaxy 
at xyz is rotated to XYZ. The origin of both coordinate systems 
corresponds to the centre of the galaxy, and xy or XY indicate 
the mid-plane of the disk of the galaxy. The Sun is located at 
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Fig. 1. Rotation of the Galaxy from xjz-coordinates to XYZ- 
coordinates. 



X = -8.5 kpc, and y = z = 0. The rotation is accomplished in 
three steps, 

- around z by ; 

- around new x (green line in Fig. 1) by ^2; 

- around Z by ^3 . 

The inclination (/) of the new Milky-Way-like galaxy is deter- 
mined by 61 and 62 as, 

cos / = - cos bo cos /o sin Oi sin 62 + cos bo sin /q cos Oi sin O2 ... 
- sinZ?o cos^2, 

where < ^1 < 2;r, and < 62 < n.ln the simulations presented 
in this paper, we set /q = 90°, bo = 0°. The inclination angle of 
the galaxy is then simplified as cos / = cos Oi cos O2. 



3. Scale-height of the synchrotron emission 

3.1. The scale-height of the Milky Way 

We simulated the synchrotron emission of the Milky Way 
Galaxy when seen edge-on at a distance of 5 Mpc. The sam- 
pling was chosen to be 03 to ensure that averaging eff'ects are 
small. The total intensity of the Milky Way Galaxy seen edge- 
on at 408 MHz is significantly diff'erent from that modelled by 
Phillipps et al. (1981, their Fig. 8) and Beuermann et al. (1985, 
their Fig. 9) on the basis of only the 408 MHz all-sky survey 
(Haslam et al. 1982). The high-latitude emission was attributed 
to a thick disk component. Phillipps et al. (1981) obtained a huge 
halo extending in z up to about 10 kpc, and Beuermann et al. 
(1985) obtained a slightly smaller value. We show an edge-on 
view of the Beuermann et al. (1985) model in comparison with 
ours in Fig. 2, where the present z-extent is about a factor of two 
smaller. This diff'erence results from the fact that the current sim- 
ulations took into account the observed local synchrotron excess 
as reviewed in Sect. 2.1. 



1 \ r 
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Fig. 2. Edge-on view of the Milky Way Galaxy at 408 MHz. The 
black contours are from the model by Beuermann et al. (1985), 
and red contours with levels of 10 K, 100 K and 300 K from 
outside to inside are from the present 3D-model. 



3.2. The scale-height obtained for inclined galaxies 

Dumke et al. (1995, 2000) proposed a scheme to derive the scale- 
height of the synchrotron emission for highly inclined galax- 
ies. The intensity profile along the major axis is projected to 
the z axis by scaling the width by a factor cos /. For an edge- 
on galaxy, the projection results in a delta- function at the centre. 
The projected profile is then smoothed by using a Gaussian beam 
corresponding to the angular resolution of the observation. The 
smoothed profile is fitted by a Gaussian to yield a new half power 
beam width (HPBW), called the eff'ective beam size. The ob- 
served disk profile of a galaxy is commonly assumed to be expo- 
nential (oc exp(-|z|/zo,e)) or Gaussian (oc exp(-z^/zl g)), which 
is convolved to the eff'ective beam size. By fitting the newly con- 
volved profile to the observation, the scale-heights (zo, e or zo, g) 
can be obtained. In this way, both the smoothing eff'ects from 
the observations and the influence of the inclined disk can be 
accounted for. This procedure was subsequently used by many 
authors. Although this method is based on reasonable assump- 
tions, there is no independent verification. 

We simulated the 4.8 GHz synchrotron emission from the 
Galaxy viewed edge-on at a distance of 5 Mpc with a sampling 
of Of 3, made stripe integrations perpendicular to the disk, and 
normalised the profile with the intensity maximum at z = 0. The 
result is shown in Fig. 3. The distribution is well fitted by an 
exponential with a scale-height of 0.74 kpc, which is close to the 
value of 0.8 kpc for the cosmic-ray electron density in our model 
(see Table 1). The intrinsic value is thus z = 0.74 kpc, which 
one expects to obtain from observations at any inclinations and 
angular resolutions. 

We simulated galaxies with inclinations from / = 90° to / = 
70° in steps of 5°, and smoothed the total intensity maps to an 
angular resolution of T, which is about that of the Eff'elsberg 
100-m telescope at 10 GHz. Stripe integrations were made in 
the same way as for the edge-on case shown in Fig. 3 and the 
results are shown in Fig. 4. 

For / = 90° the eff'ective HPBW is V as expected. For 
smaller inclination angles the HPBW gradually broadens, which 
accounts for the contribution from the disk. For / > 80°, only 
models with exponential disks yield good fits (Fig. 4). This is 
expected because the input disk of the model is exponential. For 
small inclinations, the fits with an exponential or a Gaussian disk 
are both fairly good for intensities larger than about 10% of the 
central value (Fig. 4). However, only models with exponential 
disks reproduce the intrinsic scale-heights reasonably. The scale- 
height is overestimated by about 10% at / = 80° and about 40% 
at / = 70°. The Gaussian disks always produce too large scale- 
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z (kpc) 

Fig. 4. Average intensity distribution (shown by black dots) along z for galaxies observed with a V beam. The projected and 
smoothed profiles used to obtain the eff'ective HPBW as indicated is shown in the left column. Fits from exponential and Gaussian 
disks were shown by solid and dotted lines, respectively. 




Fig. 3. Normalised intensity distribution along z for the Milky 
Way Galaxy when seen edge-on. The fitted exponential scale- 
height is zo = 0.74 kpc. 



heights. For NGC 253 with / = 78?5, Heesen et al. (2009) found 
z = 1.7 + 0.1 kpc. For this inclination, we got a scale-height 
15% larger than expected, which exceeds the quoted error, and 
indicates the inclination limit for the method. For galaxies with 
inclination angles below / = 70°, the contribution from the disk 
masks the much fainter halo emission, so that a scale-height de- 
termination is not meaningful anymore. 

For many galaxies, models with two exponential disks seem 
to give a better fit. Their typical scale-heights from A6 cm ob- 
servations are at the order of 300 pc and 1.8 kpc (e.g. Dumke 
& Krause 1998; Soida et al. 2011). This is not the case for our 



Galaxy model, where individual source complexes in the disk are 
not included. We note that in most cases the diff'erence between 
one- and two-disk models are based on intensities below about 
10% of the central values (e.g. Fig. 6 in Dumke et al. 2000). 
More galaxies observed with a high signal-to-noise ratios are re- 
quired for a thorough statistical study. 

4. Energy equipartition 

The assumption of energy equipartition or minimum energy as- 
sumption has been widely used to estimate magnetic fields in 
galaxies or supernova remnants. Beck (2005) gave a thorough 
review on this topic and a detailed discussion on its applica- 
tion. We followed Beck & Krause (2005) to calculate the average 
magnetic field with the observed synchrotron total intensity and 
spectral index (a), and the so-called K factor, which is the ratio 
between the proton and electron number density per unit energy. 
For spiral galaxies, the K factor was traditionally set to 100 (e.g. 
Niklas & Beck 1997) and more recent evidence favors this value 
(see e.g. Heesen et al. 2011). 

The true average total magnetic field in the disk, which we 

used in the model, can be represented as ^tot = -^^reg + ^ran. 
where the regular field 5reg has an exponential shape at radii 
larger than 5 kpc and a constant value inside, and the random 
field ^ran IS homogcncous with a variance of 3 yuG (Sun et al. 
2008). The radial profile of the total magnetic field calculated 
this way is shown in Fig. 5. 

We simulated the synchrotron emission for an edge-on and a 
face-on view of the Galaxy. In both cases we assessed the aver- 



4 



X. H. Sun and W. Reich: Polarisation properties of Milky- Way-like galaxies 



S 



I •» f t 



O i=90° • i=0° 



= 100 



R (kpc) 

Fig. 5. Average total magnetic field estimates based on energy 
equipartition for an edge-on (open circle) and a face-on (filled 
circle) view of the Galaxy. The K factor is indicated. The total 
magnetic field strength of the model is shown by a solid line. 



age total magnetic field based on equipartition (Beck & Krause 
2005) with a path-length of 20 kpc and 2 kpc, respectively. For 
K factor, we used the value 100, which is the currently best 
estimate. The results are shown in Fig. 5. The magnetic field 
strengths for both views are very similar, because the emission 
from the disk dominates the total intensity compared to the emis- 
sion from the halo. At a radius of 8.5 kpc, the field strength is 
consistent with the value 6+2 juG calculated by E. Berkhuijsen 
(see Beck 2001) based on the emission in the plane from the 
Beuermann et al. (1985) model (see their Fig. 6a). 

The magnetic field strength estimated from energy equipar- 
tition with a K factor of 100, however, overestimates the true 
field strength by a factor of nearly 2. To roughly reproduce 
the intrinsic total magnetic field, the K factor should be about 
10 (Fig. 5). However, as already pointed out above, ^ = 100 
seems to be well settled (Heesen et al. 201 1). The overestimated 
magnetic field strength indicates that the cosmic-ray particle en- 
ergy is about four times larger than the magnetic energy in our 
model. The energy density of the magnetic field including both 
regular and random components is eb = (B^^g + B^^JI%7i « 
0.5 X 10"^^ ergs cm"^ at = 8.5 kpc, and correspondingly 
the energy density of cosmic-ray particles is about ecre = 
2.0 X 10"^^ ergs cm~^. The local total pressure from cosmic rays 
and magnetic field is Pcre+b = I^cre ~ 1.2 x 10"^^ dyn 
cm~^. As summarised by Cox (2005), to support the gas against 
the gravity, the local non-thermal pressure from cosmic rays, 
magnetic fields and gas motions all together should be about 
2.8 x 10"^^ dyn cm~^. The dynamic pressure is therefore about 
1.6 X 10"^^ dyn cm~^, which dominates the pressure from mag- 
netic fields. Recent MHD simulations by Hill et al. (2012) show 
that the magnetic field does not play an important role in sup- 
porting the gas. 

It is clear that the equipartition magnetic field strength does 
not agree with the true value obtained directly from the model. 
This is not unexpected, as we know that our input model does 
not satisfy energy equipartition. In the same sense, it is also very 
uncertain whether spiral galaxies generally conform to the en- 
ergy equipartition. It is thus questionable whether magnetic field 
estimates based on energy equipartition are reliable. 

One may ask if it is possible to use the equipartition value for 
the magnetic field and adapt the cosmic-ray density accordingly 
to model the Galaxy. We have tested that and could not reach 
an agreement with the observations. A higher regular magnetic 
field will increase the RMs of pulsars and extragalactic sources, 
which could only be compensated by decreasing the thermal 



electron density. The latter, however, has been constrained by 
pulsar dispersion measures. Increasing the random field compo- 
nent will increase the total synchrotron emission and hence the 
cosmic-ray electron density must be reduced. This consequently 
reduces the polarisation intensity. To reproduce the observations, 
the regular magnetic field has to be increased, which is diflficult 
as explained at the beginning. 

In case the Milky Way is not exceptional, the equipartition 
magnetic field strength estimated for ^ = 100 might also be 
lower for other galaxies. The reduction by a factor of two we 
found for the Milky Way corresponds to a magnetic energy den- 
sity decrease by a factor of four and has a large influence on the 
average magnetic pressure. For example, when equipartition is 
assumed the magnetic field energy density clearly dominates the 
turbulent energy density outside of the inner kpc of NGC 6946 
(Beck 2007). This is very puzzling, because the turbulent mag- 
netic field is believed to be connected to turbulent gas motions. 
When reducing the magnetic energy density by a factor of four, 
however, it is almost equal to the turbulent energy density. 



5. Rotation measures of galaxies 

We simulated two samples of galaxies at 4.8 GHz (A6.25 cm) and 
8.4 GHz (A3.6 cm), which are standard frequencies for nearby 
galaxy observations. RM maps of nearby galaxies are mainly 
calculated from polarisation measurements at these frequencies 
(e.g. Beck 2007). It is generally assumed that the galaxies are 
Faraday thin at these frequencies, which means that the observed 
polarised emission originates from the same volume. We placed 
the simulated galaxies at a distance of 5 Mpc, corresponding to 
a maximum size of about 23' X 23', when the galaxies are seen 
face-on. The resolution or the grid size of the simulations is Of3 
or about 440 pc, which is comparable to observations and pro- 
vides some details of the large-scale emission from the galaxies. 

We simulated maps for all Stokes parameters and from them 
calculated polarisation angles (PA) at both frequencies. From 
the polarisation angle maps, we obtained RM for each individual 
pixel as. 



RM = 



PAi - PA2 



(2) 



where the subscripts stand for the two bands, e.g. A = 6.25 cm 
and A = 3.6 cm. The ^7r-ambiguity of the polarisation angles at 
these two bands corresponds to a RM of 1194 rad m~^. We as- 
sume that ^=0 and the RM values are smaller than the ambiguity, 
which is certainly true except for the centre regions. 

The Faraday depth (FD) of a galaxy in rad m~^ is defined as. 



FD = 0.81 



Jfar 



n.B\\dL 



(3) 



where Ug is the thermal electron density in cm"^ and B\\ is the 
magnetic field parallel to the line-of-sight in fiG. The integral is 
calculated along the line-of-sight from the far to the near side of 
the galaxy. It is obvious that FDs are very important to under- 
stand the properties of magnetic fields and thermal electron den- 
sities. In the case that the synchrotron-emitting medium and the 
thermal gas are uniformly mixed, one would expect FD = 2RM 
(e.g. Sokoloff' et al. 1998). However, as we will show below, the 
relation does not hold anymore when the frequency-dependent 
polarisation horizon is shorter than the size of a galaxy. 

The RM and FD maps for a nearly face-on (/ = 20°) and a 
nearly edge-on (/ = 80°) galaxy are shown in Fig. 6. In general. 
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Fig. 6. Total intensity (I) at 1.4 GHz, rotation measure (RM) calculated from 4.8 GHz and 8.4 GHz simulations, Faraday depth 
(FD), and depolarisation maps between 4.8 GHz and 8.4 GHz (DPcx) and between 1.4 GHz and 4.8 GHz (DPlc) for two Milky 
Way galaxies with inclination angles / = 20° and / = 80°, respectively. 



both RM and FD maps show the same intensity pattern, which 
means that RM maps based on observations correctly indicate 
the orientation of the large-scale field of a galaxy. The field re- 
versals in the Milky Way Galaxy and the opposite sign of the 
halo magnetic field above and below the Galactic plane are re- 
flected in the RM maps. 

The RMs multiplied by a factor of two were compared to 
the corresponding FDs in Fig. 7. For nearly face-on galaxies 
(/ = 20°) both distributions are basically consistent with each 
other, although a small amount of scatter is visible. For nearly 
edge-on galaxies (/ = 80°), however, the differences are signif- 
icant. Large deviations from FD (Fig. 7) are seen for large FD 
values originating in the disk. The reason is that the simulated 
4.8 GHz and 8.4 GHz polarised emission does not originate from 
the same volume. The polarisation horizon at 4.8 GHz is smaller 
than that at 8.4 GHz. Higher- frequency observations are needed 
to overcome this effect and provide RMs proportional to FDs. 

Structure functions (e.g. Sun & Reich 2009) for RMs and 
FDs were calculated following 



D(SO) =< [RM(0) - RM(0 -h SO)f 



(4) 



where < . . . > stands for ensemble average and S6 is the an- 
gular separation. The results are shown in Fig. 7. The structure 



functions for RM, when scaled by a factor of two, is equiva- 
lent to shifting constantly the structure functions for the original 
RM upwards along the ordinate axis in Fig. 7. For galaxies with 
high inclinations, the slopes of the structure functions differ for 
RMs and FDs, while for low-inclination galaxies RMs mimic 
FDs quite well after being multiplied by a factor of two. 

In summary, for face-on galaxies RMs derived from 4.8 GHz 
and 8.4 GHz polarisation measurements can be safely used to 
infer magnetic fields after being scaled by a factor of 2. For edge- 
on galaxies or highly inclined galaxies, however, RMs are only 
reliable in the outskirts of these galaxies, but still indicate the 
orientation of the magnetic field direction. 

In Fig. 6 we show maps of the depolarisation between the 
two simulated maps at 4.8 GHz and 8.4 GHz, which show clear 
differences related to their inclination. For the definition of de- 
polarisation (DPcx) we follow Heesen et al. (2011) as. 



DP- 



cx 



100- 



PL 



4.8 GHz 



W8.4GHZ \4.8 



8.4 



(5) 



where JS = a - 2 = -3 is the spectral index for brightness 
temperature. A value of DPcx = 100 means no depolarisation, 
and DPcx = means total depolarisation. In the depolarisation 
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Fig. 7. RMs times a factor of two versus FDs (top row). The solid line indicates the expectation for uniformly mixed thermal and 
synchrotron emission. The structure functions (bottom row) from the FD maps are shown by black dots, those from RM maps and 
RMs scaled by a factor of two by squares. 



map for the galaxy with inclination / = 20°, the values exceed 
92% everywhere. For the nearly edge-on galaxy with i = 80°, 
the DPcx along the major axis is small and the values close 
to the centre are about 40%. Although our model does not in- 
clude discrete thermal complexes such as H II regions, which are 
certainly present along the spiral arms, the depolarisation maps 
show a patchy morphology. This results from the various ran- 
dom magnetic field components along line-of- sights. The highly 
inclined galaxy NGC 253 with / = 78?5 was recently studied by 
Heesen et al. (2011) using polarisation observations at 4.8 GHz 
and 8.4 GHz. Their depolarisation map shows similar values as 
our results, although their distribution of depolarisation patches 
is more complex than ours. The main reason is the missing dis- 
crete source complexes in our model. 

We also calculated the depolarisation (DPlc) between 
1.4 GHz and 4.8 GHz (Fig. 6). For low inclination galaxies, 
DPlc is around 40% within 5' and there is a slight asym- 
metry along the minor axis direction, which is caused by the 
toroidal halo field. We note that some galaxies, for instance 
NGC 6946, show a clear asymmetric polarisation at frequencies 
around 1.4 GHz with complete depolarisation towards a large 
patch (Beck 2007; Heald et al. 2009). This might indicate a 
more complex halo field configuration than in the Milky Way. 
A quadrupole field configuration was discussed by Braun et al. 
(2010), but asymmetries could also be ascribed to some large- 
scale distortions from a regular configuration. In our Galaxy, 
discrete large RM or polarisation features exist, for instance the 



well pronounced large Fan region (e.g. WoUeben et al. 2006), 
which are commonly considered to be local. If such regions re- 
side out of the Galactic plane or in the halo, they would make the 
polarisation properties of our Galaxy asymmetric and thus simi- 
lar to others. More investigations are needed to improve Galactic 
emission models. 



6. Integrated polarisation properties of Milky Way 
type galaxies 

Spiral galaxies are believed to occupy a large fraction of the po- 
larised radio sources at 1.4 GHz at yuJy flux densities (Stil et al. 
2007). However, they are difficult to be resolved individually 
even with instruments such as ASKAP (Johnston et al. 2007) 
or the future SKA. It is therefore important to investigate the 
polarisation properties of unresolved spiral galaxies. 

We simulated 1000 Milky- Way-like galaxies by randomly 
setting 6i and O2 (see Sect. 2.2) at each of the five frequencies: 
1.4 GHz, 2.7 GHz, 4.8 GHz, 8.4 GHz, and 22 GHz. The dis- 
tance was set to 10 Mpc and the angular resolution to 18" or 
0.87 kpc. We have tried simulations at higher angular resolutions 
and the results are similar. A sample of Milky- Way-like galaxies 
at 1.4 GHz is shown in Fig. 8. 

For each object we integrated the /, U and Q maps and ob- 
tained the polarisation intensity PI = ^U^ -\- Q'^ and the polar- 
isation percentage H = PI /I. We binned the data according to 
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Fig. 8. Sample of Milky- Way-like galaxies at 1.4 GHz in total intensity. The inclination angles are marked. The position of the Solar 
System (local enhancement) is indicated by a black dot. The angular resolution is 03. The gray scale runs from to 1.5 K. 



their inclination angles with a step interval of 2°. The average 
percentage polarisation as a function of inclination is shown in 
Fig. 9. 

Face-on galaxies show a small degree of polarisation be- 
cause of the assumed symmetry of the magnetic field in the 
disk, which strongly cancels the U and Q components when in- 
tegrated. This is common for all frequencies (Fig. 9). Faraday 
eff'ects become highly significant for edge-on galaxies because 
of the large path-length along the disk emission, which renders 
the polarisation percentage low. Therefore the maximum polar- 
isation percentage is expected in between, around / = 60° at 
2.7 GHz, around / = 70° at 4.8 GHz and around / = 80° at 
8.4 GHz. The shift stems from the reduction of Faraday depo- 
larisation towards higher frequencies. At low frequencies such 
as 1.4 GHz, however, Faraday eff'ects are strong irrespective 
of inclinations, which always yields small polarisation percent- 
ages. At high frequencies, the thermal fraction in a galaxy gets 
larger, which consequently reduces the polarisation percentage 
despite reduced Faraday depolarisation and explains the lower 
percentage polarisation at 22 GHz when compared with those at 
8.4 GHz. 

The polarisation percentage at 4.8 GHz (Fig. 9) peaks at 
the same inclination range as that found by Stil et al. (2009) 
or Stil (2009), who used a geometrical galaxy model. The mod- 
els by Stil et al. (2009) all predict zero polarisation percentage 
for edge-on galaxies, which is not found in our simulations. The 
internal depolarisation for edge-on galaxies was probably over- 
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Fig. 9. Polarisation percentages of unresolved Milky Way galax- 
ies versus inclination angles at several frequencies as indicated. 



estimated by Stil et al. (2009). In those cases, the polarisation 
horizon is small (e.g. Sun et al. 201 1) and much shorter than the 
diameter of the disk. However, the disk size was used by Stil 
et al. (2009) to derive the RM fluctuations, which certainly re- 
sults in too large fluctuation and hence depolarisation. In fact, 
most of the edge-on galaxies observed at 4.8 GHz show consid- 
erable polarisation percentages (see Fig. 1 in Stil et al. 2009), 
and agree with our simulations (Fig. 9). 
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The distribution of polarisation percentages for 1000 simu- 
lated galaxies at 1.4 GHz and 4.8 GHz are presented in Fig. 10, 
which peaks at about 0.8% at 1.4 GHz and 4.2% at 4.8 GHz. 
In our distributions, the relative number of galaxies at zero po- 
larisation percentage is very small, which is clearly opposite to 
that in the histograms by Stil et al. (2009, their Fig. 6). This dif- 
ference is not caused by varying the parameters in the model of 
Stil et al. (2009), and should be explained in the following way: 
To ensure a spacial uniform distribution of galaxies, the proba- 
bility to observe a galaxy at inclination between / and / -h di is 
proportional to sin / (e.g. Stil et al. 2009), meaning the relative 
number of edge-on galaxies is maximal. Provided that edge-on 
galaxies have zero percentage as modelled by Stil et al. (2009), 
it follows that there is still a significant number of galaxies at 
zero percentage in their distributions. 

We note that the polarisation percentages can be as high as 
about 20% in the models by Stil et al. (2009) with varying mag- 
netic field properties. In this paper, we simulate Milky-Way-like 
galaxies with the same properties of the interstellar medium, as 
they were derived from observations (Sun et al. 2008), and sim- 
ply rotate the Milky Way Galaxy to realise diff'erent inclinations. 
This sample should therefore be considered as a subset of spi- 
ral galaxies. In future when more data are available from e.g. 
GALFACTS (Taylor & Salter 2010), we will try to represent 
general galaxies by modifying the Milky Way model accord- 
ingly. 
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Fig. 10. Histogram of polarisation percentages at 1.4 GHz 
(black) and 4.8 GHz (red) for 1000 simulated Milky Way galax- 
ies. 



The observed polarisation from unresolved galaxies results 
from averaging the projection of their azimuthal magnetic field 
components on the plane of the sky. Therefore we expect an 
alignment between the B -vectors (0^) and the position angles 
of the major axes (©major), which is indeed seen from the present 
simulations as shown in Fig. 1 1 . This means that polarisation an- 
gles of galaxies, when calculated from their integrated emission, 
do not vary with frequencies, implying that these galaxies do not 
show intrinsic RMs. This result has already been reported by Stil 
et al. (2009) from their simulations and is now confirmed. This 
characteristic makes unresolved spiral galaxies with optically 
measured position angles ideal background polarised sources to 
study Galactic and intergalactic large-scale magnetic fields with 
instruments like the future SKA. 




50 100 150 



Fig. 11. B -vectors (PA -h 90°) versus position angles of the ma- 
jor axes from simulated Milky Way galaxies at 1.4 GHz and 
4.8 GHz. 

7. Summary 

The Milky Way Galaxy is widely considered as a typical spiral 
galaxy. We used the observationally well constrained Galactic 
3D-emission model (Sun et al. 2008; Sun & Reich 2010) to sim- 
ulate emission of galaxies by putting the Galaxy away at any 
distances and rotating it by any angle. 

In our models, the scale-height of the synchrotron emission 
from the Milky Way is much smaller than that assumed in the 
early models (Phillipps et al. 1981; Beuermann et al. 1985). At 
that time, the evidence for the local synchrotron excess was not 
well determined and therefore not taken into account. We fit to 
our Milky Way model, when seen edge-on, an exponential scale- 
height of 0.74 kpc. 

We tested a current analysis method to derive the scale- 
height of highly inclined galaxies. We find that for Milky Way 
galaxies the scale-height apparently gets larger with decreasing 
inclination. For inclinations larger than 80°, however, the in- 
crease is below 10%. The total magnetic field estimate based 
on the equipartition assumption is about twice that of the Milky 
Way Galaxy model. 

Polarisation observations of highly inclined galaxies when 
used to derive the distribution of RMs are not reliable for their 
inner parts. At frequencies of 4.8 GHz or lower, they are not 
Faraday thin anymore. The polarised emission from two fre- 
quency observations used to calculate RMs originates from dif- 
ferent depths. This is also reflected in the slope of the structure 
functions, which clearly difl'er for RMs and FDs for galaxies 
nearly seen edge-on. 

We simulated a large sample of distant unresolved galaxies 
as it will be accessible to the future SKA and its path-finders, 
and analysed the polarisation properties of this sample for a 
large frequency range. This analysis follows that by Stil et al. 
(2009) of a large spiral galaxy sample based on their model of 
polarised emission. We qualitatively confirm the Stil et al. (2009) 
results, where the percentage polarisation drops to zero for face- 
on galaxies at all frequencies. For edge-on galaxies Stil et al. 
(2009) found almost no percentage polarisation, while we note 
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a moderate polarisation percentage decrease. The peak percent- 
age polarisation of the galaxy sample peaks at about 4.2% at 
4.8 GHz. At 1.4 GHz we find an average percentage polarisa- 
tions of about 0.8%. We confirm the previous finding by Stil 
et al. (2009) that unresolved spiral galaxies have a very well 
aligned polarisation angle in the direction of the galaxies ma- 
jor axis, which makes them ideal sources to study foreground 
Faraday rotation effects on Mpc- scales. 
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